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AND  PARAMETRIC  SOURCES  IN  SHALLOW  WATER 
T.  G.  Muir 

Applied  Research  Laboratories 
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Abstract 

Existing  theory  is  extended  to  model  the  average  field  of  sound  propa¬ 
gating  from  directive  sources  in  a  shallow  water  waveguide.  The  role  of 
narrow  beam  sources  in  the  excitation  of  modes  is  explored.  Sample  calcula¬ 
tions  are  given. 

Introduction  ’j 

'  1  2 

In  several  papers,  Weston  ’  has  advanced  a  simple  energy  flux  method, 
first  suggested  by  P.  J.  Westervelt,  for  the  purpose  of  understanding  and 
giving  perspective  to  several  characteristic  regimes  of  shallow  water  propa¬ 
gation.  Although  this  method  is  based  on  elementary  geometric  considerations, 
the  work  of  Brekhovskikh^  may  be  cited  for  precise  equivalents,  obtained  from 
first  principles,  which  support  many  of  Che  same  conclusions.  Recently, 

Shang  Er-chang^  has  used  a  different  approach  but  has  arrived  at  similar  re¬ 
sults.  All  of  these  models  are  applicable  to  linear  radiations  from  omni¬ 
directional  sources. 

The  present  paper  extends  Weston's  method  to  include  finite  beamwidth 
sources  typical  of  both  conventional  linear  and  parametric  sources  operating 
in  isothermal  water.  Although  the  results  obtained  deal  only  with  the  average 
field  intensity  as  predicted  from  consideration  of  generic  energy  flow  and 
loss  mechanisms,  they  are  nonetheless  quite  useful,  both  as  instructive  tools 
and  as  simple  models  for  engineering  design  work. 

The  incorporation  of  parametric  theory  in  the  present  paper  is  quite 
flexible  in  that  practically  any  model  that  accounts  for  parametric  generation 
in  the  freefield  provides  a  viable  starting  point  to  the  theoretical  framework. 

Linear  Source  Models 


A  shallow  water  medium  of  depth  H  overlays  and  includes  an  absorptive 
sedimentary  bottom,  having  a  critical  grazing  angle  <fic  that  defines  a  limit 
for  "total"  internal  reflection.  As  always,  <J>C  =  arc  cos  cw/cs,  where  cw  and 
cs  are  the  mean  sound  velocities  in  the  water  and  sediment  media,  respectively. 
The  water  is  assumed  to  be  shallow  in  the  acoustic  sense,  that  is,  the  ratio 
H/A,  where  A  is  the  wavelength,  is  less  than  a  few  decades. 

The  sound  source  is  aimed  horizontally  from  a  position  S  near  middepth, 
and  has  a  halfpower  beamwidth  of  2$^.  This  source  has  a  nearfield  extending 
to  a  range  R0  *»  na^/A,  where  a  is  an  effective  source  dimension.  For  the 
purpose  at  hand,  it  is  not  necessary  to  model  the  transducer  nearfield  since 
an  extrapolation  to  a  1  m  source  level  will  suffice. 
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We  begin  with  the  sonar  equation  for  the  sound  pressure  level  at  range  Rq, 
in  the  form, 

SPL  =  10  log  P  +  170.8  +  DI  -  10  log  Rq2  ,  (1) 

where  P  is  the  acoustic  power  radiated,  and  DI  is  the  directivity  index  of 
the  transducer.  Some  may  prefer  to  set  R0  =  1  m  in  these  equations;  this 
would  make  no  difference  at  long  ranges. 

At  ranges  R  >  Rq,  the  radiation  proceeds  initially  into  a  region  of 
spherical  spreading,  where  the  propagation  loss  becomes 

L  =  10  log  R2/R  2  +  a  R  ,  (2) 

°  o  w 

Here,  a  denotes  attenuation  in  the  water  column.  The  extent  of  the  spherical 
spreading  region  may  depend  on  several  conditions. 

For  fairly  large  source  beamwidths,  spherical  spreading  may  apply  only 
out  to  the  range  associated  with  the  critical  angle  of  the  bottom.  Thus,  for 
>  $c’  w^ere  is  the  halfpower  angle  of  the  beam,  we  take  the  extent  of 
this  region  as  R0  <  R  <  H/2<pc.  This  simply  means  that  energy  has  a  tendency 
to  be  radiated  into  the  bottom  at  high  grazing  angles.  The  geometry  for  large 
beam  insonif ication  is  sketched  in  Fig.  1. 


FIGURE  1.  LARGE  BEAM  INSONIFICATION 
SHOWING  SPHERICAL  AND  CYLINDRICAL  REGIMES 


For  smaller  source  beamwidths,  such  that  <pi^  <  4>c ,  spherical  spreading 
applies  out  to  the  range  at  which  the  beam  effectively  "fills"  the  water 
column,  i.e.,  over  the  span  RQ  <  R  <  H/24^.  The  geometry  for  this  case  is 
sketched  in  Fig.  2. 
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FIGURE  2.  MODERATE  BEAM  INSONIFICATION 
SHOWING  SPHERICAL  AND  CYLINDRICAL  REGIMES 

Beyond  the  region  of  spherical  spreading,  we  enter  one  of  cylindrical 
spreading,  brought  on  by  increased  bottom  reflectivity.  Provided  we  are 
above  the  cutoff  frequency  (=  cw/4H<J>c),  the  propagation  loss  becomes 

L  -  10  log[(RH/24>)/R02]+  awR  ,  (3) 

where  <J>  =  <j>c  or  fa  ,  whichever  is  smaller. 

For  large  beams,  the  cylindrical  spreading  region  is  usually  quite  limited 
in  lateral  extent  and  in  some  cases  is  insignificantly  small.  It  will  be  shown 
that  decreasing  the  source  beamwidth  increases  the  extent  of  this  region, 
normally,  this  region  is  characterized  by  the  internal  reflection  of  a  large 
number  of  rays  or  modes.  The  finite  attenuation  of  either  conceptual  entity 
upon  bottom  reflection  ultimately  leads  to  the  next  propagation  regime,  known 
as  the  mode  stripping  region. 

Although  the  source  emits  a  continuum  of  rays,  the  medium  imposes  a 
selection  process  upon  their  existence,  much  as  a  stringed  instrument  selects 
and  resonates  preferred  harmonics  of  the  notes  played.  In  the  case  of  shallow 
water  propagation  the  analogous  "harmonics"  are  the  spatially  selected  modes, 
each  of  which  consists  of  a  pair  of  upgoing  and  downgoing  waves  ricocheting 
down  the  waveguide  on  unique  angular  paths. 

Following  Weston's  approach,^"  we  assume  that  for  any  allowed  modal  or 
eigenray  angle  <J>^,  there  is  one  loss  producing  bottom  bounce  each  time  the 
horizontal  range  cycles  through  an  incremental  distance  2H/ <p .  Thus  at 
range  R,  there  are  n  =  R<p / 2H  bottom  bounces.  If  the  loss  at  each  such  event  is 
taken  as  cup  (in  units  of  dB/rad.),  then  the  bottom  loss  accumulates  as 

L  (4>,R)  =  a<J>2R/2H  =  10  log  exp  -  |  ci<J>2R/8.68H  |  .  (4) 

2 

Weston  argued  that  the  gaussian  dependence  on  p  in  the  preceding  equation, 
should  be  tied  in  some  way  to  the  critical  angle  This  was  done  by  allow¬ 

ing  the  exponent  to  take  on  the  value  tt/4  at  $  =  4>c,  which  provides  for  a 
little  more  than  3  dB  angle  dependent  bottom  loss  at  the  range  R  ”  6.8H/apc  . 
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The  propagation  loss  in  this  region  is  found  by  substituting  the  resulting 
value  for  $c  *  4>  °  V  6.8H/aR'  into  Eq.  (3),  to  yield 

L  «  10  log |r^2  a*5  /5.2R02|+  awR  .  (5) 

This  is  a  form  of  the  well  known  three  halves  law  appropriate  for  what  is 
known  as  the  mode  stripping  region.  Thus,  at  ranges  in  the  vicinity  of  6.8H/a<f>^, 
bottom  losses  have  begun  to  strip  off  the  modes  through  a  gaussian  absorption  C 
factor  which  is  dependant  on  the  steepness  of  the  eigenray  angles.  This  start¬ 
ing  range  for  the  applicability  of  Eq.  (5)  is  of  course  applicable  to  the  broad 
beam  case,  >  <J>C.  Although  it  is  difficult  to  illustrate  mode  stripping, 
a  sketch  of  a  few  modal  eigenrays  is  offered  in  Fig.  3  for  the  case  of  insoni- 
fication  in  a  fairly  narrow  angular  sector. 


FIGURE  3.  MODE  STRIPPING  FOR  THE  SMALL  BEAMWIDTH  CASE 
SHOWING  ONE  EIGENRAY  COMPLEMENT  FOR  TWO  MODES  DEPICTED 


A  similar  criterion  may  be  applied  to  the  narrow  beam  case,  tp d>c  . 
Here,  we  allow  the  exponent  in  Eq .  (4)  to  take  on  the  value  n/4  at  the  angle 
4>  «  which  provides  for  a  little  more  than  3  dB  angle  dependent  bottom  loss 
at  the  range  R  =  6 .  SH/a^2  •  The  resulting  propagation  loss  is  identical  to 
Eq .  (i ) .  The  only  difference  between  the  two  cases  is  the  range  at  which 
Eq .  (!>)  might  become  applicable. 


Since  the  stripping  of  the  higher  ordered  modes  is  accentuated  by  the 
steepness  of  their  eigenray  angles,  these  modes  are  eventually  attenuated, 
leaving  only  one  propagating  mode.  By  the  time  this  has  occurred,  the  surviv¬ 
ing  mode  will  itself  have  suffered  some  angle  dependent  attenuation.  Although 
the  sole  surviving  mode  is  sometimes  not  that  of  lowest  order  (depending  for 
example  on  bottom  peculiarities),  we  here  assume  that  it  is  the  first  mode, 
again  in  consonance  with  Weston's  analysis. ^ 


For  many  real  bottoms  with  phase  shifts  near  it  radians,  the  modal  angle 
for  the  first  mode  is  -  A/2H.  Taking  this  as  the  angle  at  which  about  half 
of  the  energy  of  the  first  mode  is  itself  stripped  off,  we  can  again  use  the 
exponent  of  Eq.  (4)  to  determine  the  range  at  which  this  occurs.  For  broad 
beams,  the  result  is  R  *  27HvaA  .  Beyond  this  range,  the  surface  and  the 
bottom  try  to  contain  a  sound  field  propagating  at  a  constant  angle,  and  we  are 
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back  to  cylindrical  spreading  plus  some  attenuation  losses.  The  adaptation 
of  Eq.  (3)  to  this  case  leads  to 

L  -  10  log  |(RH/2$1)/Ro2  )+  na^  +  awR  , 


or 

L  -  10  log  RH2/XR  2  +  RX2a/8H3  +  a  R 
o  w 


(6) 


where  the  second  term  represents  the  bottom  loss  of  the  first  mode. 

In  the  event  that  one  has  a  very  narrow  beam,  i.e.,  cpi, <  cp ^ ,  the  starting 
range  of  the  single  mode  region  is  not  determined  as  above.2  Obviously,  in 
order  to  get  such  a  highly  directive  radiation  to  propagate,  it  would  be 
necessary  to  aim  it  at'  one  of  the  eigenray  angles  appropriate  for  the  exci¬ 
tation  of  some  mode.  This  case  is  illustrated  in  Fig.  4.  In  practice,  it  is 
only  necessary  to  employ  one  beam  in  the  vicinity  of  an  eigenray  angle,  rather 
than  the  two  beams  shown.  For  the  first  mode,  this  beam  may  encompass  both 
eigenray  angles.  If  this  is  done,  it  is  reasonable  to  assume  that  modal 
propagation  would  begin  after  the  beam  has  spread  to  the  point  of  filling  up 
the  water  column.  In  this  case,  the  starting  range  becomes  R  =  H/2i(i^. 


FIGURE  4.  SINGLE  MODE  EXCITED  BY  VERY 
NARROW  BEAMS  CENTERED  ON  EIGENRAY  PAIRS 


Although  Eqs.  (2),  (3),  (5),  and  (6)  describe  successive  propagation 
regimes,  it  should  be  remarked  that  these  may  not  always  follow  in  succession. 
In  particular,  the  use  of  an  extremely  narrow  beam  appears  to  circumvent  some 
of  these  regions.  Given  that  $}j*«X/4a,  one  can  use  the  inequality  defining 
the  extent  of  the  spherical  spreading  region  to  show  that  it  would  be  skipped 
provided  a  >  2H/tt.  Similarly,  the  mode  stripping  region  would  be  skipped, 
provided 

Parametric  Source  Models 


The  application  of  parametric  arrays  to  shallow  water  (waveguide)  propa¬ 
gation  differs  in  several  respects  to  that  of  linear  sound  sources.  Physical 


Proceedings  of  The  Institute  of  Acoustics 


5.1 


Some  Simple  Propagation  Models  for  Linear  and  Parametric  Sources 
in  Shallow  Water 


and  engineering  aspects  of  this  problem  have  been  discussed  in  previous  papers'*’^ 
which  have  presented  theory  on  preferential  mode  selection  as  well  as  the 
results  of  experiments  in  both  a  laboratory  tank  and  in  a  natural  lagoon. 

Drawing  from  previous  experience,  it  is  useful  to  make  one  restriction 
concerning  the  length  of  the  parametric  interaction  volume  in  the  shallow 
water  waveguide.  Since  there  is  usually  a  phase  shift  near  it  radians  at  both 
the  surface  and  bottom  (at  low  grazing  angles),  one  must  consider  the  problem 
of  destructive  interference  of  parametric  sources  and  their  images  within  this 
volume.  This  problem  has  been  treated  for  the  freefield  case,7  and  is  illus¬ 
trated  in  this  context  in  Fig.  5.  Since  180°  phase  shifts  in  each  primary 
have  no  effect  on  the  phase  of  difference  frequency  generation,  it  can  be 
argued  that  difference  frequency  sound  generated  after  a  phase  reversal  is 
simply  180°  out  of  phase  with  that  reflected  from  the  previous  segment  of  the 
interaction  volume. 


VOLUME  VOLUME 

7 7777777777  TTTTTTTTTT7TTTTTTT7 TfTTTTTTTTT  t  /  TTTTT1  TTTTTTTm  tt'T  T)  l'TT777  Hit  t  f  fTTlfT 


FIGURE  5.  THE  PROBLEM  OF  PHASE  REVERSAL 
WITHIN  A  PARAMETRIC  INTERACTION  VOLUME 

For  the  case  of  a  shallow  water  waveguide,  it  is  convenient  to  avoid  the 
parametric  image  interference  problem  by  limiting  the  effective  length  of  the 
interaction  volume,  LA.  This  can  be  done  by  design  by  simply  restricting  LA 
to  distances  where  Rgp  denotes  the  range  at  which  the  parametric  beam 

"fills  up"  che  water  column.  For  a  source  at  middepth,  Rgp  *  H/2^.  This 
case  is  illustrated  in  Fig.  6.  Since  LA  has  an  exponential  taper,  the  definition 
of  effective  length  can  be  neither  exact  nor  critical.  Some  might  prefer 
defining  LA  as  the  1/e  absorption  distance  of  the  primaries,  i.e.,  LA  ■  l/ap, 
where  a  is  the  mean  primary  absorption  coefficient.  Although  this  is  a  safe 
and  neat  definition,  it  is  probably  overly  restrictive,  since  there  is  always 
some  bottom  bounce  loss  which  serves  to  further  truncate  the  parametric  array 
in  much  the  same  way  as  does  the  water  column  absorption.  Perhaps  a  better 
choice  is  the  -3  dB  absorption  range  LA  ■  -(in  l//2)/ap.  When  this  is  equated 
to  Rgp  one  obtains  ap5^  in  l//2)/H,  which  enables  the  primary  frequency  fp 

to  be  determined  through  ap  ■  5p(fp).  In  practice,  it  may  be  convenient  to 
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truncate  the  parametric  array  with  finite  amplitude  attenuation.  This  enables 
one  to  vary  the  array  length  by  varying  the  transmitted  power. 


FIGURE  6.  RESTRICTION  OF  PARAMETRIC  ARRAY 
LENGTHS  TO  AVOID  PHASE  REVERSAL  INTERFERENCE 


Many  theoretical  models  of  parametric  arrays  operating  in  infinite  media 
have  been  published.  The  complexity  of  these  solutions  usually  increases  with 
their  generality  and  applicability  to  various  problems.  Given  that  the  pur¬ 
pose  of  the  present  paper  is  to  address  shallow  water  applications  with  simple 
models,  we  shall  favor  the  simplest  parametric  array  descriptions. 

Most  parametric  array  models  have  their  roots  in  the  solution  to  an  in¬ 
homogeneous  wave  equation3  which  can  be  written  in  terms  of  the  difference 
frequency  sound  pressure  level  as 


Where  [•)  is  the  Green's  function,  R  is  a  vector  from  the  origin  to  the  field 
point,  r  is  a  vector  from  the  origin  to  the  element  of  integration,  V  denotes 
the  interaction  volume,  and  q  is  the  source  strength  density, 

q(r.t)  =  — §— 4  P  2  (r,t)  .  (8) 

po  co 

Here  8=1+  4(B/A),  where  B/A  is  the  parameter  of  nonlinearity,  and  pQ  and  cq 
are  the  mean  density  and  sound  speed  of  the  medium. 
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The  first,  shown  in  Fig.  7,  depicts  sound  propagating  at  500  Hz  from  a 
2.3  m  diameter  source  excited  with  80  kW  of  electrical  power,  assuming  a  50% 
transducer  efficiency.  The  water  column  is  50  m  deep,  overlying  a  sand  bottom 
having  a  sound  velocity  15%  higher  than  water,  a  critical  angle  of  30°,  and 
a  bottom  loss  a  of  12.4  dB/rad.  The  transducer  is  located  near  middepth, 
with  a  horizontal  orientation. 


RANGE  -  m 


FIGURE  7.  COMPARISON  OF  MEAN  PROPAGATION  CURVES  FOR  A 
2.3  m  diam  SOURCE  RADIATING  LINEAR  AND  PARAMETRIC  SIGNALS 
AT  500  Hz  IN  A  SHALLOW  WATER  WAVEGUIDE,  50  m  DEEP 

Since  the  linear  version  of  this  source  has  a  half  power  beamwidth  of 
80°  (ijn.  =  40°),  spherical  spreading  applies  out  to  the  critical  angle  range, 
a  distance  of  48  m.  For  this  case,  mode  stripping  becomes  significant  at  100  tn, 
while  the  single  mode  region  emerges  at  31  km.  The  3  dB  artifact  at  the 
start  of  this  region  is  discussed  in  Ref.  1.  Although  the  present  models 
treat  only  the  mean  propagation  curve,  considerable  multipath  in  the  cylin¬ 
drical  spreading  region  and  modal  interference  in  the  mode  stripping  region 
would  be  present  and  would  appear  as  wild  fluctuations  about  the  mean,  were 
they  to  be  taken  into  consideration. 

Parametric  operation  of  this  source,  also  depicted  in  Fig.  7,  arises 
from  nonlinear  interaction  of  primaries  radiated  rt  13.75  and  14.25  kHz. 

Given  the  same  available  electrical  power  and  transducer  efficiency,  one 
would  here  raJiate  about  20  kW  per  primary  tone.  The  peak  primary  source 
level  during  the  constructive  interference  phase  of  the  composite,  two  tone 
radiation  extrapolates  to  160  dB  re  1  piPa  at  1  m.  Application  of  weak  shock 
theory  indicates  this  radiation  would  form  a  discontinuity  near  the  Rayleigh 
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distance  of  40  m,  with  stronger  shocks  occurring  in  the  spherical  spreading 
field  beyond.  As  a  result,  this  parametric  array  is  truncated  by  finite  ampli¬ 
tude  attenuation  at  a  range  of  about  700  m. 

The  parametric  propagation  curve  shown  in  Fig.  7  takes  these  factors 
into  consideration,  through  appropriate  solution(s)  of  Eq.  (8).  This  includes 
a  three-dimensional  numerical  integration  in  the  primary  nearfield.  The  far- 
field  solution  was  executed  by  Professor  Rolleigh,^  who  reduced  the  problem 
to  a  single  integration  which  also  contains  a  finite  amplitude  taper  function. 
Finite  amplitude  attenuation  broadens  the  half  power  angles  of  the  parametric 
beam  to  about  ±2°,  which  leads  to  beam  intersection  with  the  bottom  at  735  m. 
Since  the  eigenray  angles  for  the  first  mode  are  centered  at  ±1.7°,  most  of 
the  energy  in  the  parametric  beam  should  couple  into  the  first  mode  of  prop¬ 
agation.  ’  The  portion  of  the  curve  beyond  the  array  truncation  and/or  beam 
intersection  point  is  therefore  constructed  from  Eq.  (7),  the  single  mode 
result  of  the  present  hiodel.  Preferential  excitation  of  the  first  mode 
under  these  conditions  minimizes  multipath  and  multimode  interference,  and 
should  make  the  mean  propagation  curve  depicted  here  a  good  representation  of 
the  actual  situation. 

Although  discrete  modal  excitation  is  an  interesting  phenomenon  that  has 
potential  as  a  research  tool  as  well  as  for  clear  channel  communications, 
there  are  situations  where  it  may  be  necessary  to  emphasize  other  aspects  of 
the  application  of  parametric  arrays  to  shallow  water  sound.  Perhaps  the 
most  important  of  these  situations  is  the  frequent  requirement  for  higher 
parametric  source  levels.  This  can  be  achieved  by  decreasing  the  ratio  of 
primary  to  secondary  frequencies.  For  the  2.3  m  transducer  in  the  50  m  wave¬ 
guide  at  hand,  we  might  not  wish  to  lower  the  primary  frequencies  as  this 
would  increase  the  length  of  the  parametric  array  which  may  lead  to  destructive 
interference  effects  within  the  interaction  volume.  A  better  approach  may 
be  to  raise  the  difference  frequency.  Choice  of  a  2  kHz  difference  frequency, 
for  example,  provides  for  a  downshift  ratio  of  7,  which  compares  to  a  value 
of  22  for  the  previous  example  at  500  Hz. 

Propagation  curves  for  a  2  kHz  example  are  shown  in  Fig.  8,  again  com¬ 
paring  linear  and  parametric  radiations.  The  half  power  angle  for  the  linear 
source  is  10°  which  is  smaller  than  the  30°  critical  angle.  This  is  in  con¬ 
tradistinction  to  the  previous  case  for  a  500  Hz  radiation.  This  diminishes 
the  vertical  angular  extent  of  steep  angle  modes  and  increases  the  range  span 
of  the  cylindrical  spreading  region.  The  mode  stripping  region  is  also  ex¬ 
panded  due  to  the  reduction  in  wavelength  and  the  consequent  excitation  of 
more  modes. 

Parametric  operation  of  the  source  again  leads  to  primary  beams  of  half 
power  angles  near  23.  Since  the  modal  angles  are  now  +n(0.43)°,  we  should 
excite  several  modes  with  the  parametric  beam.  However,  the  parametric  beam 
is  much  smaller  than  its  linear  counterpart  and  this  further  increases  tHe 
extent  of  the  cylindrical  spreading  region.  The  mode  stripping  region 
for  the  parametric  source  is  smaller  than  for  the  linear  case  simply  because 
of  the  smaller  beam  and  the  smaller  number  of  modes  available  there  for  strip¬ 
ping.  This  region  ends  at  the  same  range  for  both  linear  and  parametric  sources. 
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Within  the  cylindrical  spreading  and  mode  stripping  regions,  many  interference 
fluctuations  not  treated  by  the  present  model  will  exist.  These  should  be 
less  severe  for  parametric  excitation. 


FIGURE  8.  COMPARISON  OF  MEAN  PROPAGATION  CURVES  FOR  A 
2.3  m  diem  SOURCE  RADIATING  LINEAR  AND  PARAMETRIC  SIGNALS 
AT  2  kHz  IN  A  SHALLOW  WATER  WAVEGUIDE,  50  m  DEEP 


Discussion 


Although  this  paper  has  addressed  only  the  mean  axial  propagation  char¬ 
acter  of  linear  and  parametric  sources,  many  other  aspects  of  the  problem 
warrant  consideration.  When  the  problem  of  echo  to  reverberation  becomes 
important,  the  high  directivity  of  the  parametric  source  begins  to  weigh 
favorably  against  equivalent  broad  beam  linear  sources.  Wide  bandwidth  capa¬ 
bilities  (for  signal  processing)  as  well  as  Doppler  considerations  also  ap¬ 
pear  to  favor  parametric  sources. 

Vertical  strings  of  linear  elements  appear  to  offer  better  directivity 
than  the  linear  array  considered  here.  However,  they  are  usually  omnidirec¬ 
tional  in  the  horizontal  plane  and  are  usually  applied  to  single  mode  selection 
only  as  receivers . 

One  should  also  consider  the  electroacoustic  efficiency  of  linear  trans¬ 
mitting  arrays  in  the  very  low  frequency  regime.  In  the  preceding  discussion, 
an  efficiency  of  50%  was  assumed;  however,  this  is  usually  only  possible  with 
very  large  resonant  elements  that  WQuld  have  dimensions  of  several  meters  at 
500  Hz.  Smaller  elements  have  lower  efficiencies  and  this  may  further  favor 


Proceedings  of  The  Institute  of  Acoustics 


5.1 


Some  Simple  Propagation  Models  for  Linear  and  Parametric  Sources 
in  Shallow  Water 

the  practicality  of  the  parametric  array  for  shallow  water  applications. 

Much  work  remains  to  be  done  before  these  speculations  can  be  evaluated. 

The  present  models  are  simple  and  adequate  for  many  purposes  but  room  for 
improvement  exists.  Experimentation  is  also  in  its  infancy  and  must  of  course 
be  conducted  in  order  to  further  the  understanding  of  the  shallow  water  problem. 

Summary  of  Results 

Parameters: 

a  ■=  transducer  diameter 
X  =  wavelength 
P  =  acoustic  power 
DI  =  directivity  index 
R  =  range 

Rq  =  ita  /X  (the  Rayleigh  distance) 

SPL  =  sound  pressure  level 

pQ  ■  static  density 

k.  «*  difference  frequency  wavenumber 
a 

R  =  vector  from  the  origin  to  the  field  point 

x  =  vector  from  the  origin  to  the  element  of  integration 

V  =  volume  of  integration 

q  =  source  strength  density 

6  =  1  +  Js(B/A) 

where 

B/A  =  parameter  of ' nonlinear ity 

c  =  mean  sound  speed 
o 

t  =  time 

=  primary  pressure  field 

a  =  attenuation  coefficient  of  water 
w 

H  =  water  depth 

=  half  angle  beamwidth 
<j>  =  critical  angle  of  bottom 

a  =  bottom  loss,  dB/rad. 

L  =  propagation  loss 
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Starting  point  for  a  linear  radiation  at  a  farfield  point  R  =  7ia2/X: 

o 


SPL  -  10  log  P  +  170.8  +  DI  -  10  log  R 


(1) 


Starting  point  for  a  parametric  radiation  in  the  interaction  region: 


SPLd(R) 


(7) 


where 


q(r,u) 


4  at 


(r,t) 


(8) 


which  includes  an  appropriate  finite  amplitude  taper  function  for  the  peak 
interference  maxima  in  the  composite  primary  pressure  pp(r,t). 

Spherical  spreading  region: 


L  =  10  log  R2/R  2  +  a  R 
o  w 

valid  for  Rq  <  R  <  H/ 2<p 


Cylindrical  spreading  region: 


where  ij>  =  { 


whichever 
is  smaller 


L  *  10  log  |(RH/24>) /Rq2 j  +  c^R 
valid  for  H/2<p  <  R  <  6.8H/a<f>2  where 


Mode  stripping  region: 


whichever 
is  smaller 


L  -  10  log  [r^  H*  cx^I 5 . 2R  2  +  o  R 

o  *  w 


27H3/aX2  for  <f>^  >  ^ 


valid  for  6,8H/a$  <  R  < 


where  $ 


H/2<»», 


for  ^  Z 


c  which¬ 
ever 

*is  is 
'  smaller 
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Single  mode  region: 


L  =*  10  log  RH2/XR2  +  RX 
o 


valid  for  R  > 


3  2 

27H  /aX* 


H/2$. 


2a/8H3  +  a  R 
w 


for  «t>1 

for  <p,2!  <J> 

2  i 
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